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INTRODUCTION 


Wide  range  expressions  (interpolating  polynomials)  for  displacements  for  standard  ASTM  fracture 
testing  specimens  are  used  in  the  measurement  of  4.J-  R,  and  R  curve  determination.  They  are  used  to 
predict  crack  extension  using  the  compliance  unloading  methods.  Although  they  have  not  been  used  in 
the  measurement  of  K,c  or  for  fatigue  crack  growth  testing,  such  expressions  may  be  useful  in  these  test 
methods  and  could  be  incorporated  in  the  standards  for  these  tests  in  the  future. 

Attempts  have  been  made  in  the  past  to  develop  appropriate  expressions  (refs  1-5).  The  work  of 
Saxena  and  Hudak  (ref  1)  concentrated  on  the  compact  specimens  and  bending  samples.  The  work  of 
Kapp  et  al.  (ref  2)  attempted  to  address  all  of  the  specimens.  The  work  of  Wu  (ref  3),  Haagag  and 
Underwood  (ref  4),  and  Underwood  et  al.  (ref  5),  dealt  with  bending  samples.  The  purpose  of  this  report 
is  to  present  expressions  for  all  of  the  present  standard  specimens  that  are  similar  in  form  and  accuracy. 

PROCEDURE 


The  goals  for  the  development  of  expressions  are  to  develop  expressions  for  displacement  as  a 
function  of  crack  length  that  are  accurate  over  the  widest  possible  range  of  crack  length  and  agree  with  the 
best  available  numerical  analysis  determination  of  displacement  within  about  one  percent  or  better.  To 
develop  the  inverse  expressions  (where  the  crack  length  is  calculated  from  displacement  measurements), 
the  following  procedure  was  used.  The  developed  displacement  expression  was  used  to  determine  the 
interpolated  values  of  displacement  at  various  crack  lengths.  This  displacement  was  normalized  to  some 
form  that  can  be  used  as  a  parameter  to  predict  the  crack  length.  Three  forms  of  this  parameter  are 
presented  for  comparison:  one  suggested  by  Saxena  and  Hudak  (ref  1)  and  two  suggested  by  Kapp  et  al. 
(ref  2).  Polynomials  that  are  functions  of  the  parameters  were  then  fit  to  predict  the  relative  crack  length 
(a/W)  to  within  0.0005  W  or  better.  Such  accuracy  is  possible  for  all  of  the  specimens  considered  at  all  of 
the  displacement  locations  used  for  measurement  with  the  exception  of  the  arc  tension  sample.  This  is 
because  the  arc  tension  sample  has  many  forms  and  has  two  additional  geometric  parameters  (eccentric 
loading  ratio  (X/W)  and  radius  ratio  (ri/r2)).  For  this  sample,  it  was  not  possible  to  predict  relative  crack 
length  to  better  than  about  ±  0.003  W. 


The  nondimensional  form  for  displacement  that  will  enable  the  determination  of  accurate 
interpolating  polynomials  for  opening  displacement  as  a  function  of  crack  length  is  derived  from  the  deep 
crack  limit  and  short  crack  limits  and  was  suggested  in  Reference  2.  This  form  applies  to  crack-mouth¬ 
opening  displacements  for  three-point  bend  samples. 


E'Bvjl-a/W)2 
15.8  MW(afW) 


=  fla!W) 


(1) 


where  E'  is  the  effective  modulus  (E  for  plane  stress  and  E/(l-u2)  for  plane  strain,  (v  is  Poisson’s  ratio)),  B 
is  the  specimen’s  through  thickness,  v  is  the  total  displacement,  a  is  the  crack  length,  W  is  the  specimen 
thickness,  and  M  is  the  applied  moment  (PS/4  for  three-point  bending  samples  (where  P  is  the  applied 
load)).  This  form  has  finite  limits  as  the  relative  crack  length  (a/W)  approaches  both  the  short  crack  limit 
(a/W  =  0)  and  the  deep  crack  limit  (a/W  =  1). 

For  load-line  displacements  for  the  compact  and  disk  compact  samples  and  crack-mouth-opening 
displacements  for  arc  tension  samples,  the  form  given  by  Eq.  (1)  is  used,  with  M  replaced  by  PW  for  both 
compact  samples  and  with  P(X+W/2+a/2)  for  arc  tension  samples.  A  second  modification  for  these 
samples  is  that  the  a/W  term  in  the  denominator  is  removed.  This  term  comes  from  the  short  crack  limit 
which  is  unknown  for  the  compact  samples  and  is  too  complicated  for  the  arc  tension  sample.  For 
displacements  other  than  load-line  for  compact  specimens,  the  W  term  in  the  denominator  is  replaced  by 
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(1+A)W,  where  A  is  the  distance  from  the  load  line  to  the  point  where  the  displacement  is  being 
measured.  Finally,  for  load-line  displacements  for  three-point  bending  samples,  the  denominator  includes 
an  M2  term.  The  explicit  form  for  each  specimen  and  displacement  location  is  presented  below. 

Another  form  of  nondimensional  displacement  has  been  suggested  as  (ref  6) 

E'Bv(l-afW)2 


15. &MW(l+ajW)\a/W) 


=  fia!W) 


which  is  merely  Eq.  (1)  with  an  additional  (1+a/W)2  term  in  the  denominator.  This  term  was  added  in 
Reference  6  for  expediency  in  fitting  and  is  not  based  on  any  derived  limit  solution.  Since  the  goal  of  this 
report  was  to  determine  interpolating  polynomials  for  displacements  that  are  accurate  to  within  about  one 
percent  or  better,  the  (1+a/W)2  term  was  not  needed.  The  form  of  Eq.  (1)  was  used  for  all  of  the 
interpolating  polynomials. 

For  the  inverse  expressions,  three  parameters  have  been  suggested.  The  familiar  form  from 
Saxena  and  Hudak  is  an  empirical  parameter 


EBv  j* 


and  a/W  =  f(U).  This  parameter  has  been  used  for  many  samples.  It  has  finite  limits  as  the  value  of 
displacement  approaches  zero,  a/W  approaches  zero  and  U  approaches  one.  When  v  approaches  infinity, 
a/W  approaches  one  and  U  approaches  zero.  Although  the  parameter  has  finite  limits,  as  the 
displacement  approaches  its  two  extremes,  this  parameter  is  unlike  the  nondimensional  form  of 
displacement  of  Eq.  (1)  in  that  it  was  not  derived  from  some  known  limiting  solution.  A  derived  form  of 
an  inverse  parameter  was  presented  in  Reference  2  based  on  the  deep  crack  limit 


U  =  1  - 


15,8Af 

E'BWv 


This  form  of  U  in  Eq.  (3)  is  the  deep  crack  limit.  As  v  approaches  infinity,  U  approaches  a/W 
and  is  found  by  algebraic  manipulation  of  Eq.  (1)  without  the  a/W  term  in  the  denominator.  Although 
this  form  was  first  presented  in  Reference  2,  it  was  dismissed  as  being  an  inadequate  parameter,  since  it 
has  an  infinite  limit  as  the  displacement  approaches  zero.  The  form  chosen  in  Reference  2  as  an 
appropriate  form  was  a  combination  of  the  two  forms  presented  above 
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This  parameter  has  finite  limits  at  both  extremes.  As  v  approaches  infinity,  a/W  approaches  one 
and  U  approaches  one,  and  as  v  approaches  zero,  a/W  approaches  zero  and  U  approaches  zero.  Using 
this  form,  the  coefficients  of  the  interpolating  polynomial  are  much  simpler  than  the  coefficients  obtained 
when  using  the  form  in  Eq.  (2). 
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In  this  report,  ail  three  forms  are  used.  The  restrictions  on  the  use  of  the  parameter  in  Eq.  (3) 
(the  parameter  blows  up  when  the  displacement  and  a/W  approach  zero)  are  meaningless  since  the  forms 
of  Eq.  (2)  or  (4)  make  their  use  at  values  of  v  approaching  zero  suspect.  To  demonstrate,  the  crack- 
mouth-opening  displacements  for  the  SE(B)  sample  from  Reference  7  are  used  to  develop  values  of  the  U 
parameters  of  Eqs.  (2)  and  (3).  These  parameters  are  plotted  against  a/W  in  Figure  1.  Also  shown  in  the 
figure  are  the  third,  fourth,  and  fifth  order  least  squares  interpolating  polynomials  for  each  parameter.  It 
is  clear  that  for  both  of  these  parameters,  the  short  crack  limit  is  reached.  However,  the  manner  in  which 
the  limit  is  approached  is  inappropriate.  All  interpolating  polynomials  for  both  parameters  predict  a 
negative  crack  length  when  the  value  of  the  parameter  approaches  its  short  crack  limit.  This  means  that 
none  of  these  parameters  are  appropriate  for  predicting  the  crack  length  when  the  crack  is  small. 
Therefore,  it  was  decided  to  limit  the  applicability  of  either  of  these  two  parameters  to  predicting  the 
minimum  crack  length  for  which  there  is  accurate  numerical  determination  of  the  value  of  the 
displacement. 


Figure  1.  Comparison  of  the  U  parameters  for  SE(B)  crack-mouth-opening 
displacements  for  0  <  a/W  <  1.0. 

Since  the  two  parameters  of  Eqs.  (3)  and  (4)  were  not  used  at  very  short  crack  lengths,  the 
derived  parameter  given  as  Eq.  (3)  can  also  be  used,  since  this  parameter  only  blows  up  when  the 
displacement  is  very  small  or  when  the  crack  length  approaches  zero.  All  three  parameters  were  evaluated 
for  the  SE(B)  crack-mouth-opening  displacements  for  values  of  a/W  between  0.2  and  1.0  and  are  plotted 
in  Figure  2.  Also  plotted  in  the  figure  is  the  fifth  order  least  squares  interpolating  polynomials  for  each 
parameter.  This  plot  shows  that  the  interpolating  polynomials  are  appropriate  in  this  case. 
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Parameter 

Figure  2.  Comparison  of  the  U  parameters  for  SE(B)  crack-mouth-opening 
displacements  for  0.2  <  a/W  <  1.0. 

RESULTS  AND  DISCUSSION 


Displacements  as  a  Function  of  Crack  Lena 


The  appropriate  form  of  displacements  as  a  function  of  crack  length  is  given  in  Eq.  (1).  The 
manner  in  which  the  left-hand  side  of  Eq.  (1)  (f(a/W))  was  determined  was  to  normalize  the  numerical 
data  from  the  best  available  source  according  to  the  left-hand  side  of  Eq.  (1).  The  discreet  values  of 
f(a/W)  determined  in  this  manner  were  coupled  with  the  limit  solution  and  the  intermediate  values  of 
f(a/W)  were  determined  by  least  squares  fitting.  All  of  the  results  can  be  summarized  as  follows: 

AalW)  =  A0  *  A^ajW)  *  A2{ajW)2  +  A3(a/W)3  *  AJdW)4  *  A^dW)5  (5) 


The  explicit  form  of  f(a/W)  for  all  of  the  specimens  and  the  values  of  the  coefficients  are  given  in  Table  1. 

Some  discussion  of  these  results  is  necessary.  First,  not  all  of  these  polynomials  were  developed 
for  this  report.  In  the  table,  if  a  reference  number  appears  next  to  the  specimen  and  location  designation, 
then  the  polynomial  is  from  that  reference.  For  the  SE(B)  (CM)  displacements,  a  totally  different  form  of 
interpolating  polynomial  appears  in  E-813  and  E-1152.  This  is  the  interpolating  polynomial  from 
Reference  8.  This  polynomial  is  wrong.  It  does  not  approach  the  short  crack  limit  correctly  and  was  fit  to 
the  uncoiTected  displacement  solution  in  Reference  7.  The  corrected  displacements  for  the  SE(B)  (CM) 
were  published  in  the  errata  of  Reference  7.  Furthermore,  the  interpolating  polynomial  reported  in  the 
table  for  the  SE(B)  (CM)  was  not  determined  by  fitting  to  all  of  the  numerically  determined  displacements 
in  Reference  7.  The  displacement  at  a/W  =  0.8  in  the  errata  to  Reference  7  is  of  questionable  accuracy. 
All  of  the  other  numerically  determined  displacements  fall  smoothly  on  a  curve  that  approaches  both  the 
deep  and  short  crack  limits,  while  the  numerically  determined  value  of  displacement  at  a/W  =  0.8  clearly 
falls  off  of  the  smooth  curve.  It  was  decided  to  develop  the  polynomial  reported  in  the  table  by  excluding 
the  numerical  displacement  solution  at  a/W  =  0.8. 
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The  accuracy  column  addresses  the  accuracy  of  the  interpolating  polynomial  with  respect  to  the 
numerically  generated  displacements.  No  attempt  is  made  to  give  a  number  to  the  absolute  accuracy  of 
the  polynomials  (the  accuracy  between  the  interpolating  polynomials  and  the  true  displacement  solutions). 
The  solutions  used  for  the  fitting  are  as  follows:  SE(B)  (CM)  (ref  7);  SE(B)  (LL)  (ref  4);  C(T),  (LL),  and 
(CM)  and  D(T),  (LL)  and  (CM)  (ref  9);  C(T)  and  C(W)  (LL  +  0.1576  W)  (ref  10);  and  A(T)  (refs  2  and 
11).  It  is  clear  that  for  all  of  the  interpolating  polynomials,  the  accuracy  requirement  of  within  one 
percent  is  met  with  the  possible  exception  of  the  A(T)  samples.  As  previously  stated,  the  A(T)  sample 
includes  two  additional  geometric  variables  that  must  also  be  interpolated.  The  interpolation  in  terms  of 
eccentric  loading  (X/W)  is  not  performed.  There  are  only  two  standard  specimens  (X/W  =  0  and  X/W  = 
0.5).  To  maximize  accuracy,  a  separate  polynomial  for  each  X/W  is  included.  This  eliminated  one  of  the 
interpolating  variables.  The  other  geometric  variable  of  radius  ratio  (r2/r2)  had  to  be  addressed 
empirically.  The  simplest  way  of  accounting  for  the  rt/r2  effect  is  with  the  use  of  the  addition  of  the  linear 
function  r,/r2  as  described  in  the  table.  This  is  just  approximate,  since  the  r,/r2  effect  is  a  function  of  a/W. 
At  short  crack  lengths,  the  rt/r2  dependence  is  very  strong,  yet  at  deep  crack  lengths,  the  r,/r2  effect  is  zero. 
In  between  these  limits,  the  effect  of  r,/r2  is  nearly  constant  and  linear,  allowing  the  use  of  such  a  simple 
correction.  However,  there  is  a  price  that  is  to  be  paid,  in  that  with  the  linear  correction  reported  in  the 
table,  the  nondimensional  form  of  f(a/W)  does  not  approach  the  deep  crack  limit  for  a  specimen  with  r,/r2 
other  than  1.0.  This  means  that  the  interpolating  polynomials  for  the  A(T)  samples  are  only  accurate  over 
the  range  of  a/W  from  0.2  to  0.8. 

Crack  Length  as  a  Function  of  Displacement 

The  three  forms  for  inverting  polynomials  given  in  Table  1  have  been  used.  These  are  the 
parameters  given  as  Eqs.  (2),  (3),  and  (4).  As  with  the  polynomials  in  Table  1,  the  crack  length  is 
estimated  by  a  polynomial 

a/W  =  A0  +  AXU  +  AJJ2  +  A3U3  *  AaU 4  +  ASU5  (6) 

where  U  is  one  of  the  parameters  given  in  Eq.  (2),  (3),  or  (4).  The  major  difference  here  is  that  the 
parameter  U  is  not  based  on  the  numerical  results  themselves,  but  by  the  calculated  values  of 
displacement  given  by  the  polynomials  in  Table  1.  This  is  done  because  of  the  manner  in  which  it  is 
anticipated  that  the  expressions  will  be  used.  First,  at  a  known  value  of  a/W,  the  displacement  of  a 
specimen  will  be  measured.  This  will  enable  the  calculation  of  the  effect  modulus  E'  by  using  one  of  the 
polynomials  in  Table  1.  The  polynomials  of  the  form  of  Eq.  (6)  will  then  be  used  to  calculate  crack 
length.  Since  the  specimen  will  be  calibrated  by  the  use  of  the  interpolating  polynomials,  errors  in  crack 
length  measurement  will  be  minimized  if  the  inverse  expressions  (Eq.  (6))  are  fit  to  the  calculated  values 
of  displacement  rather  than  the  numerically  determined  values  of  displacement.  The  values  of  the 
coefficients  in  Eq.  (6)  for  each  of  the  three  forms  of  U  given  above  are  in  Tables  2,  3,  and  4. 

A  general  discussion  of  the  polynomials  are  that  the  required  accuracy  of  ±  0.0005  W  is  obtained 
using  all  forms  of  the  U  parameter  with  tb^  exception  of  the  SE(B),  load-line  displacements,  and  the  A(T) 
samples.  In  the  case  of  SE(B)  load-line  displacements,  the  total  displacement  is  the  displacement  due  to 
the  presence  of  the  crack  plus  the  deflection  with  no  crack  present.  All  of  the  U  parameters  assume  that 
the  crack  length  is  a  function  of  the  square  root  of  displacement.  For  the  SE(B)  load-line,  the  crack 
length  should  be  a  function  of  the  displacement  plus  something,  all  under  the  square  root  symbol  (ref  2). 
This  means  that  a  parameter  other  than  the  U  parameters  of  Eq.  (2),  (3),  or  (4)  should  be  used. 

However,  it  was  possible  to  predict  crack  lengths  within  0.0005  W  with  a  fifth  order  polynomial  and  U 
according  to  Eq.  (3).  Using  either  Eq.  (2)  or  (4)  for  the  U  parameter,  only  a  fourth  order  polynomial  is 
reported  in  the  tables.  This  is  because  the  use  of  a  fifth  order  polynomial  with  the  Eq.  (2)  or  (4)  form  of 
U  did  not  increase  the  accuracy,  and  a  plot  of  these  polynomials  shows  that  the  limits  are  approached  in 
an  inappropriate  manner.  For  the  A(T)  specimens,  the  errors  are  due  to  the  fact  that  the  r,/r2  effect  on 
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the  U  parameters  is  a  function  of  crack  length.  Since  we  are  attempting  to  calculate  crack  length,  we 
cannot  use  it  to  account  for  the  errors  in  U  parameter  and  then  use  the  U  parameter  to  calculate  crack 
length.  Therefore,  we  must  account  for  the  r,/r2  effect  in  an  average  way.  This  is  the  reason  for  the 
inferior  accuracy  of  the  A(T)  sample  expressions. 

Comparing  the  use  of  the  different  forms  of  U  to  predict  crack  length,  the  following  general 
observations  can  be  made.  Using  the  form  of  U  in  Eq.  (2)  results  in  inverse  polynomials  with  coefficients 
that  are  roughly  geometrically  increasing.  A*,  is  always  about  one.  A,  is  about  five.  A,  is  about  ten.  A,  is  in 
the  hundreds.  A*  is  in  the  thousands,  etc.  This  means  that  all  of  the  terms  in  the  polynomials  using  this  U 
parameter  contribute  about  the  same  to  the  calculation  of  a/W.  Accordingly,  this  indicates  that  one  must 
be  careful  when  calculating  crack  length  using  such  polynomials.  Using  the  U  parameter  of  Eq.  (4),  and 
in  particular  using  Eq.  (3),  all  of  the  coefficients  are  of  approximately  the  same  order.  In  these  cases,  the 
terms  increasing  in  order  involving  U  contribute  less  to  the  total  calculation.  For  these  cases,  the  first  few 
terms  are  more  important  than  the  latter  terms.  With  the  use  of  the  Eq.  (3)  or  (4)  forms  of  U,  most  of 
the  calculated  a/W  is  obtained  with  the  first  few  terms,  suggesting  that  these  forms  may  be  more  forgiving 
of  calculation  error. 

The  calculation  error  discussion  is  presented  because  to  calculate  a/W  using  the  inverted 
polynomials  in  Table  2,  3,  or  4,  we  do  not  use  a  directly  measured  independent  variable  to  calculate  the 
dependent  variable.  For  example,  to  use  the  interpolating  polynomials  in  Table  1,  we  measure  the  crack 
length,  divide  by  a  constant,  and  use  that  as  the  independent  variable.  When  using  the  inverse 
polynomials,  we  measure  a  slope  (v/P  or  P/v),  multiply  by  some  constants,  invert  the  result,  take  a  square 
root,  add  or  subtract  a  constant,  and  invert  the  result  just  to  get  the  independent  variable.  Many 
calculations  are  made  before  we  use  the  interpolating  polynomial  to  obtain  the  desired  result~in  this  case, 
the  crack  length.  The  possibility  for  calculation  error  in  the  evaluation  of  the  independent  variable  (the  U 
parameter)  is  great.  This  leads  me  to  argue  that  the  best  form  of  the  parameter  U  to  use  for  calculating 
a/W  is  that  form  of  U  most  resistant  to  errors  in  the  parameter  U  itself,  since  all  of  the  U  parameters  can 
be  used  to  calculate  a/W  accurately,  if  everything  is  done  properly. 

It  is  beyond  the  scope  of  this  work  to  perform  a  detailed  analysis  of  calculation  error  on  all  of  the 
interpolating  polynomials  in  Tables  2,  3,  and  4.  An  indication  of  the  results  of  such  an  error  analysis  can 
be  observed  with  a  few  numerical  experiments  on  the  use  of  all  three  forms  of  U  for  a  single  specimen. 

The  SE(B)  crack-mouth-opening  displacement  expressions  are  chosen  arbitrarily. 

It  is  assumed  that  several  types  of  calculation  errors  can  occur.  First,  it  is  observed  that  the 
coefficients  of  the  interpolating  polynomial  necessary  when  using  the  U  parameter  in  the  form  of  Eq.  (2) 
have  as  many  as  six  significant  figures.  The  other  forms  require  fewer.  In  particular,  when  using  the  Eq. 
(3)  form  of  U,  there  are  only  three  significant  figures  except  for  the  constant  term  which  has  four.  It 
should  be  pointed  out  that  any  constant  term  in  any  interpolating  polynomial  that  has  four  decimal  places, 
the  fourth  decimal  place  was  added  artificially  to  minimize  errors.  One  may  be  tempted  to  see  how 
accurate  the  expressions  are  when  all  of  the  coefficients  are  rounded  to  three  significant  figures.  Table  5 
shows  the  comparison  of  calculating  a/W  with  all  forms  of  U  and  rounded  coefficients  on  the  interpolating 
polynomials.  The  other  forms  of  calculation  error  are  rounding  of  the  calculated  U  parameter.  In  Table 
6,  the  effects  of  rounding  the  U  parameter  to  six,  four,  three,  and  two  decimal  places  are  presented.  The 
effects  of  a  relative  percentage  error  in  calculating  the  various  U  parameters  are  presented  in  Table  7. 
Finally,  the  effect  of  an  absolute  error  of  0.0001,  0.001,  or  0.01  in  the  calculation  of  U  are  presented  in 
Table  8. 


Discussing  the  effects  of  rounding  the  coefficients  first.  Table  5  presents  two  sets  of  calculations. 
The  first  data  set  uses  the  correct  form  of  U  rounded  to  four  decimal  places  (four  decimal  places  were 
used  in  the  least  squares  fitting),  and  calculates  a/W  using  interpolating  polynomials  with  all  of  the 
significant  figures  from  Table  2,  3,  or  4.  The  results  show  that  a/W  is  predicted  to  within  0.0002  W  for  the 
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Eq.  (3)  or  (4)  forms  of  U  and  within  0.0001  W  for  the  Eq.  (2)  form  of  U.  The  second  set  of  data  again 
uses  the  correctly  calculated  U  to  four  decimal  places,  but  now  with  the  coefficients  of  the  interpolating 
polynomials  rounded  to  three  significant  places.  The  errors  in  calculating  a/W  are,  of  course,  increased. 

In  this  case,  the  maximum  error  is  0.0004  W  for  the  Eq.  (3)  form,  -0.0007  W  for  the  Eq.  (2)  form,  and  - 
0.0023  W  for  the  Eq.  (4)  form.  This  means  that  the  Eq.  (3)  form  of  U  is  most  insensitive  to  the  effects  of 
significant  figures  of  the  interpolating  polynomials. 

The  effects  of  rounding  the  calculated  value  of  U  using  all  of  the  significant  figures  in  the 
coefficients  of  the  interpolating  polynomials  are  presented  in  Table  6.  There  are  three  sets  of  calculations: 
one  for  each  form  of  U.  For  each  case,  the  calculated  U  is  rounded  to  six,  four,  three,  or  two  decimal 
places.  The  first  set  of  calculations  deals  with  the  form  of  U  given  by  Eq.  (2).  The  effects  of  rounding  the 
U  parameter  are  interesting.  First,  when  rounded  to  four  decimal  places,  the  error  calculated  in  a/W  is  ± 
0.0001  W,  but  when  rounded  to  six  decimal  places,  the  error  increases  to  ±  0.0002  W.  When  rounded  to 
three  or  two  decimal  places,  the  errors  are  ±  0.0008  and  0.009  W  (maximum),  respectively.  Using  the 
form  of  U  given  by  Eq.  (4)  (the  second  set  of  data),  the  same  accuracy  ±  0.0002  W  is  obtained  with 
rounding  to  either  six  or  four  decimal  places.  The  error  is  increased  to  -0.0008  W  maximum  when  U  is 
rounded  to  three  decimal  places,  and  when  rounded  to  two  decimal  places,  the  maximum  error  is  -0.0062 
W.  When  using  the  form  of  U  given  by  Eq.  (3),  the  error  is  ±  0.0003  W  when  the  U  parameter  is 
rounded  to  either  six,  four,  or  three  decimal  places  and  is  only  -0.0028  W  maximum  when  the  parameter  is 
rounded  to  only  two  decimal  places.  These  results  show  clearly  that  the  form  of  U  parameter  that  is  least 
sensitive  to  rounding  of  the  U  parameter  is  the  form  given  in  Eq  (3). 

The  next  for'i  of  calculation  error  addressed  is  calculating  the  U  parameter  correctly  to  six 
decimal  places,  addir  either  0.1,  1,  or  5  percent,  then  rounding  to  four  decimal  places.  The  errors  in  the 
calculation  of  a/W  for  each  form  of  U  are  presented  in  Table  7.  The  first  set  of  data  uses  U  in  the  form 
of  Eq.  (2).  Miscalculating  the  U  parameter  by  0.1  percent  will  introduce  a  maximum  error  of  -0.0006  W  in 
the  calculation  of  a/W,  a  1  percent  miscalculation  results  in  an  error  of  about  -0.004  W  at  all  crack  lengths, 
and  a  5  percent  miscalculation  results  in  about  a  -0.02  W  error  at  all  crack  lengths.  With  the  form  of  U 
parameter  of  Eq.  (4),  a  0.1  percent  miscalculation  of  U  causes  an  error  that  varies  from  0.0003  to  0.001  W, 
a  1  percent  miscalculation  results  in  an  error  that  varies  from  0.005  to  0.0098  W,  and  a  5  percent 
miscalculation  results  in  an  error  that  varies  from  0.027  to  0.048  W.  Considering  the  form  of  U  given  by 
Eq.  (3),  a  relative  miscalculation  of  0.1  percent  causes  errors  that  vary  from  0  to  0.0006  W,  a  1  percent 
miscalculation  causes  errors  that  vary  from  -0.001  to  0.0049  W,  and  a  5  percent  miscalculation  causes 
errors  from  -0.005  to  0.024  W.  It  is  difficult  to  say  which  form  of  U  is  least  sensitive  to  relative 
miscalculations  of  the  U  parameter.  It  can  be  argued  that  when  using  U  in  the  form  of  Eq.  (2),  the  errors 
are  more  consistent  (they  are  roughly  equal  at  ail  crack  lengths).  For  consistency  to  be  a  concern, 
however,  it  must  be  known  that  a  miscalculation  has  been  made.  The  meter  that  must  be  used  to 
determine  which  form  is  more  resistant  to  relative  miscalculation  errors  is  to  determine  which  form  will 
give  the  ir^st  accurate  estimate  of  a/W  over  the  widest  range  of  a/W.  Using  this,  it  is  clear  that  the  form 
of  U  given  by  Eq.  (3)  is  the  most  accurate. 

The  final  kind  of  calculation  error  that  may  be  encountered  is  an  absolute  error.  Table  8 
examines  this.  The  absolute  calculation  errors  are  introduced  in  the  following  manner.  First,  the  correct 
form  of  the  parameter  is  determined,  then  added  to  it  are  either  0.0001,  0.001,  or  0.01.  The  result  is 
rounded  to  four  decimal  places,  and  the  erroneous  value  of  U  is  used  to  calculated  a/W.  Comparisons  of 
the  effects  of  these  kinds  of  errors  are  presented  in  the  same  manner  as  above.  Using  the  form  of  U  given 
by  Eq.  (2),  a  miscalculation  error  of  0.0001  causes  errors  in  the  calculation  of  a/W  of  as  much  as  -0.0005 
W,  a  miscalculation  of  0.001  in  U  causes  errors  that  vary  from  -0.0018  to  -0.0036  W,  and  a  miscalculation 
of  0.01  in  U  results  in  errors  of  between  -0.0168  and  -0.0344  W.  Using  the  form  of  U  given  by  Eq.  (4),  a 
miscalculation  of  0.0001  increases  the  accuracy  of  the  interpolating  polynomial  to  ±  0.0001  W.  A 
miscalculation  of  0.001  in  U  causes  a  uniform  error  of  about  0.0013  W,  and  a  miscalculation  of  0.01  in  U 
results  in  an  error  of  about  0.014  W.  With  the  form  of  U  in  Eq.  (3),  the  effect  of  a  miscalculation  of  U  by 
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0.0001  has  no  effect  on  the  accuracy  of  the  predicted  a/W.  A  miscalculation  of  0.001  in  the  form  of  U  of 
Eq.  (3)  causes  errors  that  vary  from  0.0003  to  0.0008  W,  and  a  miscalculation  of  0.01  results  in  errors 
ranging  from  0.0024  to  0.0075  W.  Qearly,  the  form  that  is  most  resistant  to  absolute  miscalculation  errors 
is  the  form  given  by  Eq.  (3). 

SUMMARY  AND  CONCLUSIONS 

Wide  range  expressions  have  been  developed  for  all  of  the  standard  specimens  presently  used  to 
determine  K,^  Jlc,  K-R,  and  J-R.  One  set  of  expressions  give  displacements  as  a  function  of  crack  length. 
A  second,  third,  and  fourth  group  give  crack  length  as  a  function  of  displacement  usinp  three  ways  of 
normalizing  the  displacements.  All  of  the  expressions  for  displacements  as  a  function  of  crack  length 
agree  with  numerical  displacement  solutions  to  within  about  one  percent  or  better.  Also,  all  of  the 
expressions  that  give  crack  length  as  a  function  of  displacements  will  give  calculated  values  of  a/W  within 
0.0005  W  or  better.  The  interpolating  polynomials  for  displacements  as  a  function  of  crack  length  and 
crack  length  as  a  function  of  displacement  are  used  in  tandem.  That  is,  the  inverse  expressions  for  crack 
length  as  a  function  of  displacements  were  fit  such  that  they  invert  the  interpolating  polynomials  for 
displacement  as  a  function  of  crack  length.  Since  all  three  forms  of  the  inverse  expressions  give  results  of 
about  the  same  accuracy,  an  analysis  was  performed  to  determine  which  form  of  normalized  displacement 
should  be  used  to  calculate  crack  length.  The  form  that  appears  to  give  the  best  results  is  derived  from 
the  deep  crack  limit  solution. 
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Table  2.  Interpolating  Polynomials  for  a/W  =  [(Displacement). 
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0.989  -3.463  -0.171  24.354  -72.805  84.375  10.003W,  0.2^a/W^0.8 

and  x\h2  ^  0.4 
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Table  3.  Interpolating  Polynomials  for  a/W  =  ^Displacement). 
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t  able  3.  (Continued) 


in 

(N 

ro 

o 

O 

At 

Ai 

as 


m 

r*» 

m 

VO 


OO  00 


cn 

< 


m 

vo 

oo 

vO 


Tf 

o 

d 


r» 

Ov 

fN 

rn 

vo' 

cn 

15 


LEGEND:  (CM)=  Crack  Mouth;  (IX)  =  Load  Line;  (LL+.  1 576W)  -  Load  Line  +  0. 1 576W  (V  l  location  in  E561) 


Table  4.  Interpolating  Polynomials  lot  a/W  =  {(Displacement). 
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C(W)  (LL+0.1576W)  U  =  1  -  0.302  0.423  0.282  0.217  41.429  0.205 

±0.0003W,  a/W  a  0.35 


Table  4.  (Continued) 
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Table  5.  Effects  of  Rounding  the  Coefficients  of  Interpolating  Polynomials  for  Sli(B)  Crack-Moulh-Opening  Displacement. 
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Table  6.  Effects  of  Rounding  the  Calculated  U  Parameter  on  ('rack  l  ength  Calculation  for  Sli(B)  Crack-Mouth-Opening  Displacements. 
(Columns  Headed  ’X*  Pics  are  the  U  Parameters  Rounded  to  ’X’  Decimal  Places) 
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12.36  -0.130481  0.3001  0.0001  -0.1305  0.3001  0.0001  -0.130  0.3003  0.0003  -0.13  0.3003  0.0003 

21.04  0.133492  0.4002  0.0002  0.1335  0.4002  0.0002  0.133  0.4000  0.0000  0.13  0.3987  -0.0013 

35.94  0.336934  0.5003  0.0003  0.3369  0.5002  0.0002  0.337  0.5003  0.0003  0.34  0.5019  0.0019 

64.23  0.504027  0.5998  -0.0002  0.5040  0.5998  -0.0002  0.504  0.5998  -0.0002  0.50  0.5972  -0.0028 

128.3  0.649063  0.7001  0.0001  0.6491  0.7001  0.0001  0.649  0.7001  0.0001  0.65  0.7008  0.0008 


Table  7.  Effects  of  a  Relative  Error  of  the  Calculated  U  Parameter  on  ('rack  length  Calculation  for  SE(B)  Crack-Mouth-Opening  Displacements. 
(Columns  Headed  'X'  %  are  the  Correct  U  Parameter  Plus  ’X’  Percent) 
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35.94  0.336934  0.5003  0.0003  0.3373  0.5004  0.0004  0.3403  0.5021  0.0021  0.3538  0.5096  0.0096 

64.23  0.504027  0.5998  -0.0002  0.5045  0.6002  0.0002  0.5091  0.6031  0.0031  0.5292  0.6163  0.0163 

128.3  0.649063  0.7001  0.0001  0.6497  0.7006  0.0006  0.6556  0.7049  0.0049  0.6815  0.7244  0.0244 


Table  8.  Effects  of  an  Absolute  Error  of  the  Calculated  U  Parameter  on  ('rack  l  ength  Calculation  for  SE(B)  Crack-Mouth-Opening 
Displacements.  (Columns  Headed  ’O.XXXl’  are  the  C-orrect  U  Parameter  Plus  ’O.XXXl’) 
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